High flow nasal oxygen during sedation in the cardiac catheterisation laboratory

A randomized controlled trial

Aaron Conway https://www.aaronconway.info (University of Toronto) , Peter Collins (Toronto General Hospital) , Kristina Chang (Toronto General Hospital) , Navpreet Kamboj (Toronto General Hospital) , Ana Lopez-Filici (Toronto General Hospital) , Matteo Parotto (Toronto General Hospital)

Table of Contents


Abstract

Background

High flow nasal oxygen (HFNO) permits flows of heated, humidified gas via nasal prongs at up to 70 liters per minute. Delivering oxygen at such high flow rates has multiple physiological effects that may better support the vulnerable breathing state of patients during procedural sedation. The objective of this study was to investigate the effects of HFNO in comparison to standard facemask oxygen during cardiac implantable electronic device (CIED) procedures performed with procedural sedation.

Methods A randomized controlled trial design was used with participants randomized in a 1:1 ratio to receive oxygen supplementation through a facemask as per usual practice or to receive high flow nasal oxygen during procedural sedation administered Anesthesia Assistants under the supervision of an Anesthesiologist. High flow nasal oxygen was delivered at 50 liters per minute with a 50% oxygen:air ratio. The oxygen flow rate for partcipants randomized to facemask oxygen was ≥ 10 liters per minute. Ventilation status was measured continuously using a transcutaneous carbon dioxide monitor. We pre-specified that a difference in TcCO2 >4mmHg would be clinically significant. Adverse sedation events were assessed using the Tracking and Reporting Outcomes of Procedural Sedation (TROOPS) tool. Satisfaction with sedation was measured using the Iowa Satisfaction with Anesthesia Scale (ISAS). Patient comfort with the oxygen device and Anesthesia Assistant ratings for difficulty using the oxygen supplementation device and maintaining oxygenation were measured using ordinal scales. Outcomes were analysed using Bayesian statistical models.

Results

The difference in peak TcCO2 between groups was estimated to be 0.0 mmHg (-1.3, 1.4). The estimated effect did not exceed the clinical significance threshold in either direction with high probability and there was no discernable trend observed in how the effect varied with procedure time. Particpants rated their satisfaction with sedation similarly between groups (-0.1 (-0.3, 0.2)). The probability that patients are more likely to rate comfort with the oxygen supplementation device higher with HFNO compared to the facemask is 68%. The odds of Anesthesia Assistants being more likely to rate the use of the oxygen supplementation device as more easy to use in the facemask group was estimated to be 3 times that of the HFNO group (0.3 (0.1, 0.8)). The odds of Anesthesia Assistants being more likely to rate their ability to maintain oxygen saturations as easy in the facemask group was 10 times that of the HFNO group (0.1 (0.0, 0.3)). Minor adverse sedation events related to airway and breathing were more frequent in the HFNO group.

Conclusion

There is no clear advantage for using HFNO in preference to facemask oxygen during CIED procedures performed with sedation.

Introduction

Cardiac implantable electronic device (CIED) procedures are commonly performed with procedural sedation.(Conway et al. 2014) Prophylactic oxygen supplementation is usually administered to reduce risk of hypoxemia arising from sedation-induced hypoventilation and apnea.(A. Conway, Collins, et al. 2019; Conway et al. 2013) High flow nasal oxygen (HFNO) is a promising device for oxygen supplementation in critical care and anesthesia.(Drake 2018) HFNO devices allow for heated, humidified gas with a titratable oxygen:air ratio to be administered via nasal prongs at up to 70 liters per minute. Delivering oxygen supplementation at such high flow rates has multiple physiological effects that may better support the vulnerable breathing state of patients during procedural sedation compared with the standard approach, which is typically a facemask or nasal cannula. In particular, one of the proposed physiological effects of HFNO is that it facilitates active gas exchange during times of apnea due to the highly turbulent supraglottic flow vortices.(Hermez et al. 2019) The effects of the potential disadvantages of using HFNO during sedation should also be evaluated. It is possible that the potential gains arising from delivering oxygen supplementation through the HFNO device may be offset by reduced ability to monitor ventilation from capnography waveforms when it is being used due to exhaled carbon dioxide concentrations being “washed out” by the high flows of gas. Furthermore, recent guidelines from the American Society of Anesthesiology have stated that there is insufficient evidence regarding which methods of supplemental oxygen administration (e.g., nasal cannula, face mask, or specialized devices such as HFNO) are more effective.(Apfelbaum et al. 2018) The objective of this study was to investigate the effects of HFNO in comparison to standard facemask oxygen during cardiac implantable electronic device procedures performed with procedural sedation.

Methods

Design

A randomized controlled design was used with participants randomized in a 1:1 ratio to the following treatment conditions:

Participants

Adults undergoing an elective cardiac implantable electronic device procedure with sedation administered by an Anesthesia Assistant (de novo and replacement/revision procedures) were included. As per the hospital’s policy for conscious sedation administered by Anesthesia Assistants, these patients were determined by the Anesthesiologist to be appropriate to receive care from the Anesthesia Assistants.

Exclusion criteria

  1. Under 16 years of age.
  2. Underlying condition requiring chronic oxygen supplementation.
  3. Diagnosed respiratory condition with confirmed current hypercapnia defined as confirmed PaCO2 during current admission over 45mmHg.
  4. Pre-existing untreated pneumothorax.
  5. Transesophageal echocardiography planned for the procedure.
  6. Active nasal bleeding.
  7. Complete nasal obstruction.
  8. Recent upper airway surgery or base of skull fracture.
  9. Previous participation in the study.

Interventions

All randomized participants received usual care in regard to the medications used for sedation and physiological monitoring or other interventions to support respiratory function that are considered necessary to be initiated during the procedure by the clinicians. As per the hospital policy for sedation administered by Anesthesia Assistants, oxygen therapy is administered as indicated and prescribed to maintain oxygen saturation greater than 93%. In this study, Anesthesia Assistants followed this policy. Also in accordance with this hospital policy, the Anesthesia Assistants administer a combination of sedation to target the level of conscious sedation, which is defined as “A drug-induced depression of consciousness during which patients respond purposefully to verbal commands, either alone or accompanied by light tactile stimulation.” As such, by practicing in accordance with this hospital policy, the level of sedation being administered by the Anesthesia Assistant for participants included in this trial was ‘standardized’. The hospital policy does not stipulate any dose ranges of the different classes of sedative medications to be used and there were no restrictions on the type or dose of sedation used by Anesthesia Assistants as part of this trial. The actual doses of sedation were recorded. Only the device used for oxygen supplementation delivery differed between groups in the following ways:

Facemask oxygen supplementation

Supplemental oxygen through a facemask with the flow rate chosen by the clinician responsible for sedation as per their standard practice.

High flow nasal oxygen

The Optiflow device (Fisher and Paykel Healthcare, Auckland, New Zealand), heated breathing tube and chamber, and nasal cannula was used. This system is a humidifier with an integrated flow generator, able to humidify respiratory gases and deliver them down a heated breathing tube and through the nasal cannula interface. The gas temperature was set to the ‘High’ setting (ranges 30-32º Celsius) and titrated downwards if the patient complained of irritation. The gas flow rate was commenced at 30 liters per minute prior to sedation administration and titrated up to 50 liters per minute as tolerated by the patient after sedative medication was administered. The fraction of oxygen in the gas was commenced at 50% but could be titrated according to patient requirements (i.e. increased if there is evidence of hypoventilation, airway obstruction or inadequate oxygenation, decreased during use of diathermy). Anesthesia Assistants at the site were provided with training in the use of this mode of oxygen delivery prior to study commencement by the research team.

Concomittant care

There were no restrictions on concomittant care. For examle, Anesthesia Assistants were permitted to use whichever devices for physiological monitoring that they deemed were required and titrate sedation according to their usual practice. For ventilation monitoring, Anesthesia Assistants elected to use capnography regardless of whether supplemental oxygen was delivered via HFNO or a standard facemask. The standard facemask had an integrated CO2 sampling line that could be connected to the the sidestream capnography sampling port of the anesthesia machine. For participants randomized to HFNO, Anesthesia Assistants were able to use the integrated CO2 sampling adapter in the latest model of the HFNO nasal cannula for the majority of patients (all those recruited after September 2019). Prior to this newest model being available, Anesthesia Assistants used capnography to monitor ventilation by placing a facemask with an integrated CO2 sampling line over the HFNO nasal cannula. Oxygen supplementation was delivered through the HFNO nasal cannula and CO2 was sampled from the sampling line integrated into the facemask.

Outcomes

The selection of outcomes was informed by recommendations from the Sedation Consortium on Endpoints and Procedures for Treatment, Education and Research (SCEPTER) to assess differences between groups regarding the safety and efficacy of sedation.(Williams et al. 2017; Ward et al. 2018)

The primary outcome was peak transcutaneous carbon dioxide ( TcCO2) concentration. Secondary outcomes were:

  1. Mean transcutaneous carbon dioxide concentration measured throughout the whole procedure using the VSign 2 sensor.
  2. Trajectory of transcutaneous carbon dioxide concentration as a function of time throughout the whole procedure using the VSign 2 sensor.
  3. Area under the curve of oxygen desaturation: This is a composite measure comprising the incidence, depth, and duration of oxygen desaturation events. It is calculated as the difference between the threshold (90%) and actual oxygen saturation (SpO2) summed every second during which oxygen saturation was below threshold).
  4. Adverse sedation events: measured using the Tracking and reporting outcomes of procedural sedation (TROOPS) tool at the end of procedures.
  5. Patient satisfaction with sedation: The Iowa Satisfaction with Anesthesia Scale (ISAS) will be used to measure patient satisfaction with the anesthetic used during the procedure performed.
  6. Patient comfort of oxygen delivery: Participants will be asked to rate at the end of procedures their perceived overall comfort with the oxygen delivery device used during the procedure using a 6-level ratings scale.
  7. Anesthesia Assistant rating of difficulty maintaining oxygenation status.
  8. Anesthesia Assistant rating of difficulty using oxygen delivery device.

Data collection

Instruments

Transcutaneous carbon dioxide concentrations were measured continuously using the Sentec Digital Monitoring system with VSign 2 sensor. TcCO2 monitoring provides continuous, accurate (mean bias -0.1 mmHg) and precise (95% limits of agreement within 6 mmHg) estimates of arterial CO2 (PaCO2) when the sensor is placed on the earlobe.(A. Conway, Tipton, et al. 2019) TcCO2 monitoring may provide even more precise estimates of changes in PaCO2 (mean bias 0.03 mmHg, 95% limits of agreement -0.44 to 0.38 mmHg).(Fingleton et al. 2017) The RA attached the Sentec VSign 2 sensor to the participant’s forehead using a multi-site attachment ring. Once the Sentec Digital Monitoring System displayed a stabilized TcCO2 level, the monitor was covered with a drape so that it was not visible to research staff or clinicians. The monitor was not used by the clinicians to guide treatment. Data were downloaded after procedures for analysis.

Percentage of hemoglobin saturated with oxygen (SpO2) was measured continuously using the Sentec Digital Monitoring system with VSign 2 sensor. We filtered SpO2 measurements if a pulse rate was not detectable.

Adverse sedation events were measured using the tracking and reporting outcomes of procedural sedation (TROOPS) tool.(Roback et al. 2018) Completion of the tool requires identification and description of the adverse event, the intervention, the outcome and the overall severity of the adverse event. The Anaesthesia Assistant completed this tool at the end of procedures.

Satisfaction with sedation was measured using the Iowa Satisfaction with Anesthesia Scale (ISAS). The ISAS is a questionnaire that can be used to measure patient satisfaction with the anesthetic used during a procedure performed with sedation.(Dexter, Aker, and Wright 1997; Dexter and Candiotti 2011) It contains 11 items and takes 4-5 minutes to complete. Comfort associated with oxygen delivery: Participants were asked to rate at the end of procedures their perceived overall comfort with the oxygen delivery device used during the procedure using a 6-level rating scale with ratings of ‘maximal discomfort’, ‘very uncomfortable’, ‘uncomfortable’, ‘comfortable’, ‘very comfortable’ and ‘maximal comfort’.

Adverse effects of delivering gas at a high flow rate through the nasal passages was assessed. The RA inspected intervention group participants’ skin integrity around the nasal region at the end of procedures to assess for these potentially anticipated as well as unanticipated adverse effects.

Anaesthesia Assistants were asked to rate their: 1) perceived level of difficulty in maintaining oxygenation using a 6-level rating scale with ratings of ‘extremely difficult’, ‘very difficult’, ‘difficult’, ‘easy’, ‘very easy’, ‘extremely easy’; and 2) perceived level of difficulty using the oxygen delivery device: The Anaesthesia Assistant will be asked to rate their perceived level of difficulty using the oxygen delivery device using a 6-level rating scale with ratings of ‘extremely difficult’, ‘very difficult’, ‘difficult’, ‘easy’, ‘very easy’, ‘extremely easy’.

Supplemental oxygen use and flow settings were documented by Anesthesia Assistants in a case report form provided to them by the Research Assistant. Information on the FiO2 setting, flow rate and temperature were recorded in the HFNO group. The oxygen flow rate was recorded in the control group. The Anesthesia Assistant also recorded the total doses of sedative medications used during procedures.

Sample size calculation

We estimated based on data from our prior work(A. Conway, Collins, et al. 2019) that the peak TcCO2 level in the control group will be 47 mmHg and standard deviation in both groups will be 7 mmHg. Assuming a type I error rate of 5%, a sample of 130 participants would achieve 90% power to detect a reduction in mean TcCO2 levels of 4 mmHg in the intervention period. A difference in TcCO2 levels of 4 mmHg was selected for this sample size calculation because it was used to power previous randomized controlled trials of the effect of oxygen supplementation on ventilation status in other populations, with the authors noting that an effect of this magnitude was of physiological significance.(Perrin et al. 2011; Wijesinghe et al. 2012) Also, differences in CO2 levels of a similar magnitude have been detected in previous trials evaluating the efficacy of interventions to improve sedation safety.(Baulig et al. 2015; Carmi et al. 2011; De Oliveira Jr et al. 2010, 2014; Smith et al. 2015)

Statistical analyses

All outcomes were analysed using Bayesian statistical models. Results from Bayesian models including 95% credible intervals are far easier to interpret than frequentist confidence intervals and p-values. Data and code required to reproduce the analyses is available to access in the online supplementary information. A detailed summary of the statistical models, including their prior specifications, is presented in the Appendix. Prior distributions were chosen to be weakly informative in the absence of information concerning the likely values of model parameters (Gelman et al. 2008).

The primary outcome, peak TcCO2, was compared between groups using a robust regression model. Covariate adjustments for the stratification variables obstructive sleep apnea (OSA) status and whteher or not the procedure was a cardiac resynchronization therapy (CRT) device implant were made in all models. Baseline TcCO2 concentration is included as a covariate in the TcCO2 models, as is recommended for this type of design(Egbewale, Lewis, and Sim 2014), and it is modelled using splines. A subgroup analysis was performed by expanding the model to include treatment interactions with the covariates OSA and CRT. The results of the subgroup analysis were compared to those from the primary analysis with only main effects to investigate the robustness of the conclusions to model specification.

The secondary outcome mean TcCO2 concentration was analysed using the same model used to analyze peak TcCO2. A functional analysis of variance (ANOVA) model was used to further investigate how mean TcCO2 concentration levels differ between groups as a function of procedure time.(Yue et al. 2019) The remaining continuous outcomes, average ISAS score and SpO2 AUC, were analyzed using the same model without a covariate adjustment for Baseline TcCO2.

A logistic regression model was used for adverse sedation events measured using the TROOPS tool. A proportional-odds model was used for ordinal outcomes including participant ratings of comfort with the supplemental oxygen device, and Anesthesia Assistant ratings for difficulty maintaining oxygenation status and their rating of difficulty using supplemental oxygen device. The secondary outcome of mean TcCO2 was compared between the intervention and control groups using a Bayesian one-way functional analysis of variance (ANOVA) model .

Posterior inference for all models except the functional ANOVA model was performed using Hamiltonian Monte Carlo through the brms package (Bürkner 2017), version 2.12.0. For this set of models, 2000 posterior samples were obtained from 4 independent chains of 2000 samples, where the first 1000 warm-up samples were discarded. Posterior inference for the functional ANOVA model was performed using the Integrated Nested Laplacian Approximation (Rue, Martino, and Chopin 2009) through the INLA package, version 20.4.18. The marginal posterior distribution of parameters were summarized by their mean and a 95% credible interval defined by the interval spanning the 2.5% and 97.5% percentiles of their distributions. The clinical significance of treatment effects relating to TcCO2 concentration were evaluated by computing the posterior probability that an effect exceeds 4 mmHg in either direction.

In the case of missing outcome data, we first assumed the data to be missing completely at random (MCAR) [cite Rubin] and inference was performed using a complete-case analysis. When the proportion of missing data was large and the MCAR assumption was unlikely to be satisfied, a sensitivity analysis was performed to investigate the robustness of the conclusions of the complete-case analysis.

Results

Participants

A CONSORT flow diagram is presented in Figure 1. From August 2019 to March 2020, we screened 270 patients undergoing CIED procedures. A total of 130 participants were randomized. Although we initially planned (as outlined in the protocol) to recruit up to 150 patients, the decision was made to complete recruitment due to the onset of clinical research restrictions at the institution in response to COVID-19. One participant was subsequently excluded from the study because their procedure was cancelled after the randomization was performed. One further participant, who was randomized to the HFNO group, had their procedure rescheduled to a time that the Research Assistant was not available. As such, this participant received oxygen via standard face mask and TcCO2 data were not collected. For two participants, the TcCO2 sensor failed to callibrate prior to commencement of the procedure, so they were not able to be included in the analyses for the TcCO2 outcomes. Fourteen of the Anesthesia Assistants reported that it was the first time they had used HFNO. Most (n=29; 45%) reported having used HFNO between 2 and 5 times with only one reporting using the device more than 10 times.

CONSORT flow diagram

Figure 1: CONSORT flow diagram

A summary of demographic and clinical characteristics of the participants is presented in Table 1. The sample was mostly elder and male. Anesthesia Assistants’ rated the ASA Physical Classification Status for participants as either III or IV, reflecting the underlying cardiovascular disease and multiple comorbidities. Obstructive sleep apnea was common, with 27% of participants reporting a diagnosis of this condition. About 20% of procedures were for cardiac resnchronisation therapy.

Table 1: Participant Characteristics
Characteristic High Flow nasal oxygen, N = 641 Face mask oxygen, N = 651
Age (years) 67 (14) 70 (13)
Gender
Female 19 (30%) 17 (26%)
Male 45 (70%) 47 (72%)
Prefer not to say 0 (0%) 1 (1.5%)
Other 0 (0%) 0 (0%)
Smoking history
Never 23 (36%) 25 (38%)
Current 7 (11%) 7 (11%)
Past 34 (53%) 33 (51%)
Obstructive sleep apnea 17 (27%) 18 (28%)
Uses CPAP 9 (14%) 12 (18%)
Admission source
Ward 17 (27%) 18 (28%)
Day surgery 45 (70%) 44 (68%)
CVICU 2 (3.1%) 3 (4.6%)
CICU 0 (0%) 0 (0%)
ASA classification status
I 0 (0%) 0 (0%)
II 0 (0%) 0 (0%)
III 21 (33%) 16 (25%)
IV 43 (67%) 49 (75%)
Procedure
PPM 11 (17%) 17 (26%)
PPM generator change 6 (9.4%) 7 (11%)
PPM lead revision 0 (0%) 1 (1.5%)
ICD 19 (30%) 13 (20%)
ICD generator change 10 (16%) 13 (20%)
ICD lead revision 2 (3.1%) 0 (0%)
CRT-D 11 (17%) 11 (17%)
CRT-P 2 (3.1%) 2 (3.1%)
Wound revision 0 (0%) 0 (0%)
Other 3 (4.7%) 1 (1.5%)
Charlson Comorbidity Index 4.46 (2.14) 5.15 (2.51)
Total dose of midazolam (mg) 1.58 (0.84) 1.45 (0.71)
Total dose of propofol (mg) 100 (126) 88 (104)
Total dose of fentanyl (mcg) 71 (28) 76 (56)

1 Statistics presented: mean (SD); n (%)

Comparisons between groups

Primary outcome

Results of comparisons between HFNO and facemask oxygen are presented in Table 2. TcCO2 concentrations for all patients throughout procedures are displayed in Figure 2, with the longest procedure highlighted as a reference. The effect of HFNO on the peak TcCO2 was estimated to be 0.0 mmHg (-1.3, 1.4) and the probability that it exceeds the 4 mmHg clinical significance threshold in either direction is extremely low. The estimates for the treatment effect in the baseline, OSA, and CRT subgroups when treatment interactions were included in the model did not exceed the clinical significance threshold in either direction with high probability. Results for subgroup analyses are available in the online supplementary information.

Table 2: Results of intention-to-treat analyses
Outcome Summary value Randomization Effect type Estimated treatment effect (95% CI)
Face mask oxygen High flow nasal oxygen

Peak TcCO2

N

65

61

Mean (sd)

49.0 mmHg (6.9)

47.8 mmHg (9.7)

Absolute difference

0.0 mmHg (-1.3, 1.4)

Mean TcCO2

N

65

61

Mean (sd)

44.3 mmHg (5.9)

42.7 mmHg (7.2)

Absolute difference

-0.1 mmHg (-1.3, 1.1)

SpO2

N

65

63

SpO2 <90% event

6 (9%)

7 (11%)

Median (IQR) Area under SpO2 desaturation curve

50 (10, 226.5)

32 (4, 56)

ISAS score

N

63

63

Mean (sd)

2.1 (0.9)

2.0 (1.0)

Absolute difference

-0.1 (-0.3, 0.2)

Patient comfort

N

65

63

Maximal comfort

17

9

Odds ratio

1.2 (0.6, 2.2)

Very comfortable

13

26

Comfortable

22

22

Uncomfortable

10

3

Very uncomfortable

2

2

Maximal discomfort

1

1

Difficulty maintaining oxygenation status

N

31

52

Extremely easy

17

9

Odds ratio

0.1 (0.0, 0.3)

Very easy

10

14

Easy

4

17

Difficult

0

6

Very difficult

0

4

Extremely difficult

0

2

Difficulty using oxygen delivery device

N

32

52

Extremely easy

17

15

Odds ratio

0.3 (0.1, 0.8)

Very easy

9

17

Easy

6

20

Minor airway or breathing event

N

65

64

Yes

2

9

Odds ratio

6.4 (1.3, 43.0)

No

63

55

Odd ratios are interpreted as the odds of the event occuring in the HFNO group compared with the odds of the event occuring in the facemask group.
TcCO2 = Transcutaneous carbon dioxide concentration
SpO2 = Percentage of hemoglobin saturate with oxygen
ISAS = Iowa Satisfaction with Anesthesia Scale
95% CI = 95% credible intervals

Secondary outcomes

The effect of high flow nasal oxygen on the mean TcCO2 concentration of the whole procedure was estimated to be -0.1 mmHg (-1.3, 1.1) and the probability that it exceeds the 4 mmHg clinical significance threshold is again low. Differences in TcCO2 concentration level between groups over time is presented in Figure 2. The estimated effect did not exceed the 4 mmHg clinical significance threshold in either direction with high probability and there is no discernable trend observed in how the effect varies with procedure time. Precision decreases as procedure time increases, reflecting the shrinking number of participants to compare at those times.

Transcutaneous carbon dioxide measurements throughout procedures

Figure 2: Transcutaneous carbon dioxide measurements throughout procedures

The odds ratio for a minor adverse airway or breathing outcome, as measured by the TROOPS tool, for the HFNO group compared with the facemask group was estimated to be 6.4 (1.3, 43.0). This effect estimate is very imprecise due to the small number of events.

The effect of HFNO on average ISAS score was estimated to be -0.1 (-0.3, 0.2), where the effect is measured in units of absolute difference. Average ISAS scores from two participants were not calculated due to missing responses for at least one of the component ISAS items.

The odds ratio for Anesthesia Assistant ratings of difficulty maintaining oxygenation status and difficulty using the oxygen delivery device as estimated using a complete-case analysis are 0.1 (0.0, 0.3) and 0.3 (0.1, 0.8), where a value less than 1 indicates a greater level of difficulty for respondents in the HFNO group. It may be simpleter for interpretation to reverse the terms. For example, the odds of Anesthesia Assistants being more likely to rate the use of the oxygen supplementation device as more easy to use (i.e., “extremely” or “very easy” versus “easy”) in the facemask group was estimated to be 3 times [i.e., 1/0.33] that of the HFNO group (0.3 (0.1, 0.8)). The odds of Anesthesia Assistants being more likely to rate their ability to maintain oxygen saturations as easy in the facemask group was 10 times [i.e., 1/0.1] that of the HFNO group (0.1 (0.0, 0.3)). It should be noted, however, that the Anesthesia Assistant ratings of difficulty using the oxygen device and difficulty maintaining oxygenation were missing 45 and 46 responses, respectively, likely due to the survey being voluntary. It is unlikely that missingness among these ratings occurred completely at random, so a best- and worst-case imputation approach was used to investigate the impact that the missing data could have on the results in extreme cases. The best- and worst-case sensitivity analysis gave estimates ranging between 0.0 (0.0, 0.1) and 3.3 (1.7, 6.6) for difficulty maintaining oxygenation status and from 0.1 (0.0, 0.2) and 5.0 (2.5, 10.4) for difficulty using the oxygen delivery device. These estimates suggest the directionality of the effect could be positive or negative with high probability, so conclusions of the complete-case analysis are not robust to assumptions about the values for the missing data.

A visualization of the amount of missing data within SpO2 trajectories is available in the online supplementary information. Due to the considerable amount of missing data observed in the SpO2 data, comparisons of the SpO2 AUC between groups was not performed.

Oxygen flow rates

Participants randomized to the HFNO group received flow rates at 50 litres per minute or higher for the majority of the time (Figure 3). Two participants who were randomized to HFNO did not receive this intervention at all during procedures at the discretion of the Anesthesia Assistant, with the rationale that the high flow of oxygen interfered with capnography monitoring. Four participants who were randomized to HFNO stopped receiving this intervention at a certain timepoint during procedures at the discretion of the Anesthesia Assistant, again with the rationale that capnography monitoring was not sufficient with the HFNO device.

Oxygen flow rate in litres per minute for each participant

Figure 3: Oxygen flow rate in litres per minute for each participant

Discussion

The primary result of this trial is that HFNO delivered at flow rates of 50 litres per minute to patients undergoing elective CIED procedures with sedation is highly unlikely to decrease or increase peak TcCO2 concentration by a clinically important amount, which we considered would be a difference of more than 4mmHg. As such, although a prior physiological modelling study identified a mechanism by which HFNO promotes carbon dioxide clearance,(Hermez et al. 2019) the magnitude of any such effect was insufficent to produce an important difference in ventilation status during procedural sedation. This result is consistent with prior clinical research in the non-sedation context. The difference in PaCO2 observed between HFNO (5.81 kPa; sd=1.1) and standard facemask oxygen (5.6 kPa; sd=1.0) from a randomized trial of 20 patients who were receiving pre-oxygenation for induction of anesthesia prior to emergency surgery was not significant (p=0.631).(Mir et al. 2017) Likewise, in a larger trial of pre-oxygenation with 80 patients, the end-tidal CO2 in the first breath after intubation was not significantly different between HFNO (5.0 kPa; sd=0.8) and standard facemask (5.3 kPa sd=1.0) oxygen supplementation (p=0.18).(Lodenius et al. 2018) Importantly, in contrast to these trials where ventilation status was assessed at one specific point in time with either PaCO2 or ETCO2 samples, we used continuous TcCO2 monitoring so that we could estimate differences in ventilation between groups over the whole duration of procedures. There was no discernable trend observed in how the effect varied over time, lending strength to the assertion that ventilation status is not impacted by the use of HFNO. Based on the evidence at hand regarding the impact of HFNO on ventilation, we can infer that using this device during procedural sedation is not likely to prevent potential adverse sedation events that may be precipitated by hypercarbia. In addition, authors of the physiological modelling study(Hermez et al. 2019) proposed that an improvement in oxygen transport would presumably go hand‐in‐hand with the enhancement of carbon dioxide clearance. The lack of difference in TcCO2 between groups in our study would suggest that the mechanism of any protective effects against hypoxemia during sedation that may be observed with the use of HFNO is therefore more likely to be due to superior pre-oxygenation(Lodenius et al. 2018; Patel and Nouraei 2015) rather than active gas exchange during hypoventilation.

Capnography is widely considered to be an essential aspect of physiological monitoring during sedation.[Conway, Douglas, and Sutherland (2016); conway2019pre; conway2014clinical] One concern about using HFNO during sedation is that the ability to monitor capnography waveforms is diminished and, as a consequence, sedated patients may have more prolonged, undetected episodes of hypoventilation when it is used.(Horiuchi and Graham 2020) If this were the case, presumably we would have observed higher TcCO2 concentrations in the HFNO group. We did not observe this to be the case in our study. Therefore, the use of HFNO during sedation should not be precluded solely on the basis that capnogrpahy monitoring is ineffective in this context.

The evidence base for the effects of HFNO therapy for procedural sedation is currently limited. One large(Lin et al. 2019) and three small randomized controlled trials of HFNO during procedural sedation have been published in 2019, with several more on-going trials registered on clinical trial registries.(Eugene et al. 2020) The primary outcomes for all the trials to date have focused on investigating the imapact of HFNO on oxygenation. Results have been inconsistent. One of the small trials randomized 60 participants undergoing bronchoscopy to receive HFNO at 50 liters per minute with 100% oxygen or to receive oxygen at 10-15 liters per minute through a facemask.(Douglas et al. 2018) There was no difference observed between the treatment groups for the primary outcome, which was the proportion of patients who experienced oxygen desaturation (defined as SpO2 90%). Another trial randomized 59 morbidly obese patients undergoing endoscopy to receive a fraction of inspired oxygen concentration of 0.36 either via HFNO at a flow rate of 60 liters per minute or via nasal cannula at 4 liters per minute.(Riccio et al. 2019) Again, there was no difference in the primary outcome of oxygen desaturation (SpO2 90%). The third study randomized 30 participants undergoing dental sedation into three groups to receive a fraction of inspired oxygen concentration of 0.4 either via HFNO at a flow rate of 50 liters per minute, via HFNO at a flow rate of 30 liters per minute or via nasal cannula at 5 liters per minute.(Sago et al. 2015) Participants randomized to the HFNO groups had higher nadir blood oxygen levels recorded than the low flow oxygen group. In contrast, a large trial of 1994 participants undergoing gastroscopy with propofol sedation reported a large reduction in risk of hypoxemia (8.4% in the control group and 0% in the HFNO group).(Lin et al. 2019) This resuslt is likely explained by the large difference in FiO2 that was delivered between the two groups. In the HFNO group participants received 60 liters of 100% oxygen per minute and in the control group participants received just 2 liters of oxygen per minute.

Results of our study have revealed some interesting insights into the impact that using HFNO during sedation for CIED procedures has on patient-reported outcomes. Overall patient satisfaction with sedation is very likely to be similar regardless of the type of oxygen supplementation device used, with no difference observed between groups in ISAS scores. There is a reasonably high probability (68%) that patients are more likely to rate comfort with the oxygen supplementation device higher with HFNO compared to the facemask. As such, there does not seem to be any negative impact on patient-reported outcomes associated with using HFNO during sedation.

In contrast to these mostly postive findings from the patient perspective, we identified that the HFNO device was rated as more difficult for Anesthesia Assistants to use compared with the standard facemask. However, it should be noted that none of the Anesthesia Assistants rated the HFNO device as difficult to use and most had very limited experience using the device. Most participants reported they had used HFNO between 2 and 5 times. Potentially, increased experience with HFNO would influence clinicians’ perceptions about the difficulty using the device in comparison to standard facemask oxygen.

There is also a high probability (84%) that minor adverse sedation events related to airway and breathing are more likely to occur when HFNO is used. The suspected etiology noted for these events was oxygen desaturation. It is possible that the default combination for the flow rate (50 liters/minute) and oxygen:air blend (50% O2) settings in the HFNO group was not equivalent to the amount of oxygen supplementation received in the facemask group. Most participants in the facemask group received >10 liters per minute of 100% O2. Therefore, to reduce the number of minor adverse sedation events related to airway and breathing for patients receiving HFNO, clinicians could consider using a higher setting for the oxygen:air blender.

Participants in our study mostly received combinations of propofol, midazolam and fentanyl for sedation, which is a common and recommended regimen for CIED procedures.(Furniss and Sneyd 2015) There is the potential that different results may be observed for patients who receive higher doses of sedation or who experience more frequent or prolonged episodes of hypoventilation and apnea than those in our study. It should also be noted that participants randomized to HFNO received oxygen at flows at or exceeding 50 liters per minure for the majority of the time. HFNO devices can deliver higher flow rates. It is of course unknown if results would have differed considerably if higher flow rates were used.

Oxygen desaturation is not a common event when oxygen supplementation at flow rates between 6-10 liters per minute through a face mask during procedures performed with sedation in the cardiac catheterisation laboratory.(Conway et al. 2013, 2014) Results from trials of HFNO during sedation conducted in other clinical settings, such as bronchoscopy and gastrointestinal endoscopy, where desaturation occurs more often, can therefore not be directly generalised to the cardiology context.

Limitations

Ventilation status was measured using TcCO2 monitoring, which is not a perfect substitute for PaCO2. The primary outcome was peak TcCO2 and we accounted for the correlation between baseline and peak measurements by including the baseline measurements as a covariate in the model. However, an inherent problem with this approach is that it is potentially sensitive to how the baseline and peak measurements were chosen. We chose the first second of TcCO2 at the time sedation was first administered as the baseline measurement and the maximimum TcCO2 observed over the whole procedure as the peak. It is unknown if results are robust to a different baseline measurement, such as the average of a certain number of seconds around this timepoint.

A further potential limitation is that we did not blind participants or clinicians to group assignment. We included only patients undergoing elective CIED procedures. Extrapolation of our findings to other procedures performed in procedural settings or emergency procedures where the patient characteristics may be very different to our sample is not recommended. The small dropout and cross-over rate is unlikely to have exerted a major impact on our estimates of the effect of HFNO during CIED procedures performed with sedation. Based on the data collected about the oxygen device settings used in both HFNO (flow rate, oxygen:air blend) and facemask (flow rate) groups, we do not suspect there was any significant bias associated with suboptimal application of the intervention in either group.

Conclusion

We identified no clear advantage for using HFNO in place of facemask oxygen during CIED procedures performed with sedation. Ventilation, as measured by TcCO2 monitoring, did not differ by a clinically important amount between groups. There is a higher probability that patients will be more comfortable during procedures with HFNO in comparison to the facemask, but overall patient satisfaction with sedation was similar between groups. Anesthesia Assistants rated the HFNO device as more difficult to use than facemask oxygen.

Appendix

Robust regression

A robust regression model is used for the analysis of continuous outcomes peak TcCO2, mean TcCO2, average SPO2, and average ISAS score. Let \(y_i\) be the value of the continuous response for the \(i\)th patient, \(x_{Baseline,i}\) be a baseline covariate, and \(x_{HFNO,i}\) be an indicator variables which is nonzero if the patient belonged to the HFNO treatment group. Associated with each patient are additional stratification indicator variables for OSA status \(x_{OSA,i}\) and CRT status \(x_{CRT,i}\). We supress the patient index \(i\) on patient covariates to simplify the notation.

The likelihood of the robust regression model is \[\begin{equation} y_i | \nu, \mu_i, \sigma \sim \text{Student-t}\left(\nu, \mu_i, \sigma\right), \end{equation}\] where \(y_i\) is the response value of the \(i\)th patient, \(\text{Student-t}\left(\nu, \mu, \sigma\right)\) is a generalized Student’s t-distribution with degrees of freedom \(\nu\), location parameter \(\mu\) and scale parameter \(\sigma\). The location parameter for the \(i\)th patient is assumed to depend on their covariates as

\[\begin{equation} \mu_i = f(x_{\text{Baseline}}) + \beta_0 + \beta_1 x_{HFNO} + \beta_2x_{OSA} + \beta_3x_{CRT}, (\#eq:robust_link) \end{equation}\]

where \(f\) is a non-linear continuous function of a continuous covariate \(x_{Baseline}\) and \(\beta_j, j=0,1,2,3\), are regression coeficients. For peak and mean TcCO2, \(x_{Baseline}\) is taken as the first recorded value of TcCO2 of the patient during their procedure. For average SPO\(_2\) and average ISAS score there is no corresponding baseline covariate and so \(f\) is omiited from @ref(eq:robust_link). Prior distributions of model parameters are chosen to be diffuse in the absence of information about their likely values. In the following list of prior specifications, \(\text{Normal}(\mu,\sigma)\) denotes a normal distribution with mean \(\mu\) and standard deviation \(\sigma\). * \(f\): The function are estimated using a thin-plate spline basis and the random-effect formulation used in brms [cite]. The spline basis dimension is 20 and the standard deviation parameter governing the random-effect distribution is given a \(\text{Student-t}\left(3, 0, 10\right)\) prior. * \(\beta_0\): intercept coefficient is given a \(\text{Normal}(\mu,25)\) prior, where the location parameter \(\mu\) is set equal to the mean value of the response. * \(\beta_1, \beta_2, \beta_3\): coefficients of indicator variables are given a \(\text{Normal}(0,25)\) prior. * \(\nu\): degrees of freedom parameter is given a \(\text{Gamma(2,0.1)}\) prior proposed for such parameters in the absence of any information (cite). * \(\sigma\): scale parameter is given a \(\text{Student-t}\left(3, 0, 10\right)\) prior.

Logistic and proportional odds regression

A logistic regression model is used for the adverse affect “patient experienced minor respiratory event”. A proportional-odds model with logit link is used for ordinal outcomes patient comfort of oxygen delivery, and Anesthesia Assistant rating of difficulty maintaining oxygenation status and rating of difficulty using oxygen delivery device. The link function in each model is written as a function of covariates in the same way as in the robust regression case (). The coeficients \(\beta_0, ... \beta_3\) are given a \(\text{Normal}(0,10)\) prior distributions in the absence of information about their likely values.

Functional analysis of variance

The functional observations \(\mathbf{y}_{ij}\) are taken to be the vector of TcCO2 measurements of the \(i\)th patient belonging to randomization group \(j\) measured at 1 second intervals from the procedure’s start until the procedure’s end. The group variable \(j=1\) indicates membership to the face mask oxygen group and \(j=2\) indicates membership to the high flow nasal oxygen group. Associated with each patient are stratification indicator variables for their OSA status \(x_{OSA}\) and CRT status \(x_{CRT}\). All functional observations are aligned to have time \(t=0\) occur at the start of each procedure. The time resolution of observations is reduced to measurements every 30 seconds to reduce the computational burden of model fitting.

The one-way functional ANOVA model assumes that the TcCO2 measurments \(\mathbf{y}_{ij}\) for the \(i\)th patient has the likelihood function \[\begin{equation} y_{ij} | \boldsymbol{\eta}_{ij},\Sigma_i \sim \text{Normal}\left(\boldsymbol{\eta}_{ij}, \Sigma_i\right), \end{equation}\] where \(\boldsymbol \eta_{ij}\) is the mean vector and \(\Sigma_i\) is a covariance matrix. The mean vector is assumed to have the form \[\begin{equation} \boldsymbol\eta_{ij} = \boldsymbol\mu + \boldsymbol\alpha_j + \beta_1x_{OSA} + \beta_2x_{CRT} + b_i, \end{equation}\] where \(\boldsymbol\mu\) is the baseline functional effect common to all patients, \(\boldsymbol\alpha_j\) is the \(j\)th treatment group functional effect,\(\beta_1\) and \(\beta_2\) are effects for OSA status and CRT status, respectively. A normally-distributed random intercept \(b_i \sim N(0,\sigma_b^2)\) is included to account for different average levels in TcCO2 observed between patients. The baseline constraint parameterization of \(\alpha_1 = 0\) is used to ensure that the model is identifiable. Under this parameterization, \(\boldsymbol\mu\) is the mean level of TcCO2 of the face mask oxygen group and \(\boldsymbol\mu + \boldsymbol\alpha_2\) is the mean level of TcCO2 of the high flow nasal oxygen group.

The set of TcCO2 measurements from a patient form a time-series that is inadequately modelled by assuming independent and identically distributed normal errors in (3). Mis-specification of the covariance structure can lead to credible intervals for the treatment effect that are artificially narrow. To better model the error structure of the data, the covariance matrix \(\Sigma_i\) is assumed to have the form of an autoregressive process of order 1 with the same set of correlation and variance parameters for every patient, denoted \(\rho\) and \(\sigma^2_{AR1}\), respectively.

The model is scaled to have generalized variance of 1 to aid in prior specification (Sørbye and Rue 2014). The prior distrbutions of model parameters are chosen to be diffuse in the absence of information about their likely values. They are as follows:

Apfelbaum, J, JB Gross, RT Connis, M Agarkar, DE Arnold, CJ Coté, and A Tung. 2018. “Practice Guidelines for Moderate Procedural Sedation and Analgesia 2018: A Report by the American Society of Anesthesiologists Task Force on Moderate Procedural Sedation and Analgesia, the American Association of Oral and Maxillofacial Surgeons, American College of Radiology, American Dental Association, American Society of Dentist Anesthesiologists, and Society of Interventional Radiology.” Anesthesiology 128 (3): 437–79.

Baulig, Werner, Marija Keselj, Barbara Baulig, Sandra Guzzella, Alain Borgeat, and Jos’e Aguirre. 2015. “Transcutaneous Continuous Carbon Dioxide Tension Monitoring Reduced Incidence, Degree and Duration of Hypercapnia During Combined Regional Anaesthesia and Monitored Anaesthesia Care in Shoulder Surgery Patients.” Journal of Clinical Monitoring and Computing 29 (4): 499–507.

Bürkner, Paul-Christian. 2017. “Brms: An R Package for Bayesian Multilevel Models Using Stan.” Journal of Statistical Software 80 (1): 1–28. https://doi.org/10.18637/jss.v080.i01.

Carmi, Uri, Mordechai R Kramer, Dmitry Zemtzov, Dror Rosengarten, and Oren Fruchter. 2011. “Propofol Safety in Bronchoscopy: Prospective Randomized Trial Using Transcutaneous Carbon Dioxide Tension Monitoring.” Respiration 82 (6): 515–21.

Conway, Aaron, Peter Collins, Kristina Chang, Sebastian Mafeld, Joanna Sutherland, and James Fingleton. 2019. “Sequence Analysis of Capnography Waveform Abnormalities During Nurse-Administered Procedural Sedation and Analgesia in the Cardiac Catheterization Laboratory.” Scientific Reports 9 (1). https://doi.org/10.1038/s41598-019-46751-2.

Conway, Aaron, Clint Douglas, and JR Sutherland. 2016. “A Systematic Review of Capnography for Sedation.” Anaesthesia 71 (4). Wiley Online Library: 450–54.

Conway, Aaron, Karen Page, John Rolley, and Paul Fulbrook. 2013. “Risk Factors for Impaired Respiratory Function During Nurse-Administered Procedural Sedation and Analgesia in the Cardiac Catheterisation Laboratory: A Matched Case–Control Study.” European Journal of Cardiovascular Nursing 12 (4): 393–99.

Conway, Aaron, John Rolley, Karen Page, and Paul Fulbrook. 2014. “Trends in Nurse-Administered Procedural Sedation and Analgesia Across Cardiac Catheterisation Laboratories in Australia and New Zealand: Results of an Electronic Survey.” Australian Critical Care 27 (1). Elsevier: 4–10.

Conway, Aaron, Elizabeth Tipton, Wei-Hong Liu, Zachary Conway, Kathleen Soalheira, Joanna Sutherland, and James Fingleton. 2019. “Accuracy and Precision of Transcutaneous Carbon Dioxide Monitoring: A Systematic Review and Meta-Analysis.” Thorax 74 (2): 157–63.

De Oliveira Jr, Gildasio S, Paul C Fitzgerald, Nora Hansen, Shireen Ahmad, and Robert J McCarthy. 2014. “The Effect of Ketamine on Hypoventilation During Deep Sedation with Midazolam and Propofol: A Randomised, Double-Blind, Placebo-Controlled Trial.” European Journal of Anaesthesiology 31 (12): 654–62.

De Oliveira Jr, G S, S Ahmad, P C Fitzgerald, and R J McCarthy. 2010. “Detection of Hypoventilation During Deep Sedation in Patients Undergoing Ambulatory Gynaecological Hysteroscopy: A Comparison Between Transcutaneous and Nasal End-Tidal Carbon Dioxide Measurements.” British Journal of Anaesthesia 104 (6): 774–78.

Dexter, Franklin, John Aker, and Will A Wright. 1997. “Development of a Measure of Patient Satisfaction with Monitored Anesthesia Care the Iowa Satisfaction with Anesthesia Scale.” Anesthesiology 87 (4): 865–73.

Dexter, Franklin, and Keith A Candiotti. 2011. “Multicenter Assessment of the Iowa Satisfaction with Anesthesia Scale, an Instrument That Measures Patient Satisfaction with Monitored Anesthesia Care.” Anesthesia & Analgesia 113 (2): 364–68.

Douglas, N, I Ng, F Nazeem, K Lee, P Mezzavia, R Krieser, D Steinfort, L Irving, and R Segal. 2018. “A Randomised Controlled Trial Comparing High-Flow Nasal Oxygen with Standard Management for Conscious Sedation During Bronchoscopy.” Anaesthesia 73 (2): 169–76.

Drake, Matthew G. 2018. “High-Flow Nasal Cannula Oxygen in Adults: An Evidence-Based Assessment.” Annals of the American Thoracic Society 15 (2): 145–55.

Egbewale, Bolaji E, Martyn Lewis, and Julius Sim. 2014. “Bias, Precision and Statistical Power of Analysis of Covariance in the Analysis of Randomized Trials with Baseline Imbalance: A Simulation Study.” BMC Medical Research Methodology 14 (1): 49.

Eugene, Axelle, Lucie Fromont, Adrien Auvet, Olivier Baert, Willy-Serge Mfam, Francis Remerand, Thierry Boulain, and Mai-Anh Nay. 2020. “High-Flow Nasal Oxygenation Versus Standard Oxygenation for Gastrointestinal Endoscopy with Sedation. The Prospective Multicentre Randomised Controlled Odephi Study Protocol.” BMJ Open 10 (2). British Medical Journal Publishing Group.

Fingleton, James, Steve Mckinstry, Janine Pilcher, Mark Weatherall, Richard Beasley, and George Bardsley. 2017. “Accuracy of Transcutaneous Carbon Dioxide Measurement for Change over Time” 22: 98.

Furniss, Stephen S, and J Robert Sneyd. 2015. “Safe Sedation in Modern Cardiological Practice.” Heart 101 (19). BMJ Publishing Group Ltd; British Cardiovascular Society: 1526–30.

Gelman, Andrew, Aleks Jakulin, Maria Grazia Pittau, Yu-Sung Su, and others. 2008. “A Weakly Informative Default Prior Distribution for Logistic and Other Regression Models.” The Annals of Applied Statistics 2 (4). Institute of Mathematical Statistics: 1360–83.

Hermez, L. A., C. J. Spence, M. J. Payton, S. A. R. Nouraei, A. Patel, and T. H. Barnes. 2019. “A Physiological Study to Determine the Mechanism of Carbon Dioxide Clearance During Apnoea When Using Transnasal Humidified Rapid Insufflation Ventilatory Exchange (Thrive).” Anaesthesia 74 (4). Wiley: 441–49. https://doi.org/10.1111/anae.14541.

Horiuchi, Akira, and David Y Graham. 2020. “High-Flow Nasal Cannula Oxygen Therapy Obscures Severe Respiration Depression, Providing an Illusion of Safety.” Gastrointestinal Endoscopy 91 (3). Elsevier: 720–21.

Lin, Yuxuan, Xiaoqing Zhang, Lizhi Li, Mengyun Wei, Bin Zhao, Xiaojing Wang, Zhiying Pan, Jie Tian, Weifeng Yu, and Diansan Su. 2019. “High-Flow Nasal Cannula Oxygen Therapy and Hypoxia During Gastroscopy with Propofol Sedation: A Randomized Multicenter Clinical Trial.” Gastrointestinal Endoscopy 90 (4). Elsevier: 591–601.

Lodenius, Å, J Piehl, A Östlund, J Ullman, and M Jonsson Fagerlund. 2018. “Transnasal Humidified Rapid-Insufflation Ventilatory Exchange (Thrive) Vs. Facemask Breathing Pre-Oxygenation for Rapid Sequence Induction in Adults: A Prospective Randomised Non-Blinded Clinical Trial.” Anaesthesia 73 (5). Wiley Online Library: 564–71.

Mir, F, A Patel, R Iqbal, M Cecconi, and SAR Nouraei. 2017. “A Randomised Controlled Trial Comparing Transnasal Humidified Rapid Insufflation Ventilatory Exchange (Thrive) Pre-Oxygenation with Facemask Pre-Oxygenation in Patients Undergoing Rapid Sequence Induction of Anaesthesia.” Anaesthesia 72 (4). Wiley Online Library: 439–43.

Patel, A, and SAR Nouraei. 2015. “Transnasal Humidified Rapid-Insufflation Ventilatory Exchange (Thrive): A Physiological Method of Increasing Apnoea Time in Patients with Difficult Airways.” Anaesthesia 70 (3): 323–29.

Perrin, Kyle, Meme Wijesinghe, Bridget Healy, Kirsten Wadsworth, Richard Bowditch, Susan Bibby, Tanya Baker, Mark Weatherall, and Richard Beasley. 2011. “Randomised Controlled Trial of High Concentration Versus Titrated Oxygen Therapy in Severe Exacerbations of Asthma.” Thorax 66 (11): 937–41.

Riccio, Christina A, Stephen Sarmiento, Abu Minhajuddin, Dawood Nasir, and Amanda A Fox. 2019. “High-Flow Versus Standard Nasal Cannula in Morbidly Obese Patients During Colonoscopy: A Prospective, Randomized Clinical Trial.” Journal of Clinical Anesthesia 54: 19–24.

Roback, M G, S M Green, G Andolfatto, P L Leroy, and K P Mason. 2018. “Tracking and Reporting Outcomes of Procedural Sedation (Troops): Standardized Quality Improvement and Research Tools from the International Committee for the Advancement of Procedural Sedation.” British Journal of Anaesthesia 120 (1): 164–72.

Rue, H., S. Martino, and N. Chopin. 2009. “Approximate Bayesian Inference for Latent Gaussian Models by Using Integrated Nested Laplace Approximations.” Journal of the Royal Statistical Society 71 (2): 319–92. http://myaccess.library.utoronto.ca/login?url=https://search-proquest-com.myaccess.library.utoronto.ca/docview/37127083?accountid=14771.

Sago, Teppei, Nozomu Harano, Yuki Chogyoji, Masahito Nunomaki, Shunji Shiiba, and Seiji Watanabe. 2015. “A Nasal High-Flow System Prevents Hypoxia in Dental Patients Under Intravenous Sedation.” Journal of Oral and Maxillofacial Surgery 73 (6): 1058–64.

Smith, Stephen B, Shawn Carr, Stacey Psikula, Anita Das, and Katherine Grichnik. 2015. “A Pilot Study on the Effect of Nasal Continuous Positive Airway Pressure on Arterial Partial Pressure of Carbon Dioxide During Spinal Anesthesia with Intravenous Sedation for Total Knee Arthroplasty.” Anesthesia & Analgesia 120 (2): 479–83.

Sørbye, Sigrunn Holbek, and Håvard Rue. 2014. “Scaling Intrinsic Gaussian Markov Random Field Priors in Spatial Modelling.” Spatial Statistics 8: 39–51. https://doi.org/https://doi.org/10.1016/j.spasta.2013.06.004.

Ward, Denham S., Mark R. Williams, John W. Berkenbosch, Maala Bhatt, Douglas Carlson, Phillip Chappell, Randall M. Clark, et al. 2018. “Evaluating Patient-Centered Outcomes in Clinical Trials of Procedural Sedation, Part 2 Safety: Sedation Consortium on Endpoints and Procedures for Treatment, Education, and Research Recommendations.” Anesthesia & Analgesia 127 (5). Ovid Technologies (Wolters Kluwer Health): 1146–54. https://doi.org/10.1213/ane.0000000000003409.

Wijesinghe, Meme, Kyle Perrin, Bridget Healy, Mark Weatherall, and Richard Beasley. 2012. “Randomized Controlled Trial of High Concentration Oxygen in Suspected Community-Acquired Pneumonia.” Journal of the Royal Society of Medicine 105 (5): 208–16.

Williams, Mark R., Denham S. Ward, Douglas Carlson, Joseph Cravero, Franklin Dexter, Jenifer R. Lightdale, Keira P. Mason, et al. 2017. “Evaluating Patient-Centered Outcomes in Clinical Trials of Procedural Sedation, Part 1 Efficacy: Sedation Consortium on Endpoints and Procedures for Treatment, Education, and Research Recommendations.” Anesthesia & Analgesia 124 (3). Ovid Technologies (Wolters Kluwer Health): 821–30. https://doi.org/10.1213/ane.0000000000001566.

Yue, Yu (Ryan), David Bolin, Håvard Rue, and Xiao-Feng Wang. 2019. “Bayesian Generalized Two-Way Anova for Functional Data Using Inla.” Statistica Sinica 29 (2): 741–67.